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lonospheric interactions 
with EME signals 


By Giorgio IK1 UWL and Flavio IK3XTV 


Date: 2012-aug-03 — Station IK1UWL - Band 144 MHz 


QRG_ DF DT dB] UTC 

144.143 -129 0-2 0 1.7| 63 4 -22}0008 CQ OX3LX HP15 

144.143 -138 3-1 0 1.7] 64 3 -20]0010 CQ OX3LX HP15 

144.143 -144 0-2-1 1.7] 48 3 -22}0012 CQ OX3LX HP15 

144.143 -153 0-2-1 1.9) 43 5 -22)/0014 CQ OX3LX HP15 

144.143 -161 1-1 0 1.5} 31 1 -21]0016 CQ OX8LX HP15 

144.143 -170 0-1 0 1.7) 15 2 -21]0018 F6EHVK OX3LX HP15 OOO 
144.143 -176 0 0 0 3.6) 0 4-18/0020 RRR 

144.143 -185 1 0 0 1.7/18 5 -20]0022 RK3FG OX3LX HP15 OOO 
144.143 -191 0 0 0 1.0) 0 4-19/0026 RRR 

144.143 -199 1-1 0 1.7] 10 4-19]0028 CQ OX3LX HP15 

144.143 -205 0-2 0 1.5] 4 4-18]0030 I3MEK OX3LX HP15 OOO 
144.143 -214 0 0 0 3.6) 0 3-18/0032 RRR 

144.143 -217-1 0-1 2.1] O 4-18)0034 IZ3KGJ OX3LX HP15 OOO 
144.143 -226 0 0 0 1.0) 0 4-23/0038 RRR 

144.143 -229 -1 -1-1 1.6 0040 CQ OX8LX HP15 

144.143 -232 -2 -1-1 1.8 0042 CQ OX8LX HP15 

144.143 -238 0-1-1 1.8 0044 CQ OX3LX HP15 

144.143 -243 3-1 1 1.8 3 -20]0046 IK1UWL OX3LX HP15 OOO 
144.143 -246 0 0 0 1.0 4 -18}0048 RRR 


MAP65 can be a research tool. 
Besides decodes for ham activity, 
it measures also levels and polarity. 


With this tool we started to research what 
happens in the ionosphere. 
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signal level (dB) 


Level oscillated from -18 to -22 dB 


- Partially ionized gas layer between ~50 and ~1000 km 
height and permeated by Earth’s magnetic field is a 


turbulent ocean, roughened by high speed winds. 
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: = Jet stream in the troposphere, always from west 
The winds « Jet stream in the mesosphere (reverses seasonally) 


=" zonal currents in the stratosphere with seasonal inversion 


1 
o° 70° 90° 
: Equator 
Summer hemisphere Winter hemisphere 


Courtesy: Istituto di Scienze dell'Atmosfera e del Clima — CNR, Italia. 
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LU. lonospheric effects: Attenuation, Deviation, Rotation of the wave 
LJ Winds cause undulations and waves (TIDs), so free electron density 


varies in space and time. 


LI These fluctuations of electron density have a lens effect on our signals, 


focusing or defocusing them. 

- Moon is wide 0.5 degrees 

- our beam is wide many degrees 
- change of width changes gain 


The Travelling lonospheric Disturbances 
(TIDs) 


Cl Horizontal — Horizontal 
ass wavelenght phase velocities 


LSTIDs | 51000 Km 0.5..3h | 300..1000 mis 
Large scale 


MSTIDs | 499.4000 Km |12 min...1h | 100..300 mis 
Medium scale 


SSTIDs <100 Km A few <200 m/s 
Small Scale minutes 


Source: INGV Istituto Nazionale di Geofisica e Vulcanologia - Italy 


1500 | \ \ Fs 9,9) MHz/ / 


1400 
1300 
1200 
1100} 
1000 PX. 
900 
800 & 
7005 
600 
500 
400 
300 


Defocusing 


Distance (km) 


6) 


2400 2200 2000 1800 1600 1400 1200 1000 800 600 400 200 O 
North—South distance (km) 


36 38 40 42 48 50 52 34 


44 46 
Latitude (°) 


Image source:Research and Technology Organisation. North Atlantic Treaty Organisation. 
Characterising the lonosphere. Author: G. Wyman (January 2009) 


LI Fast scintillations caused 
by lunar libration 
and ionospheric turbulence 


(ssTIDs, periods of minutes) 


LY Slower fluctuations from msTIDs 
(observed at mid latitudes every day) 
(300 km wavelength, wind 100 m/s 


= 360 km/h, period 50 minutes) 
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Band dependence (ionospheric 


refraction is proportional to 1/f? ) 


Regions dominating the effect 
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Fig. 8. Detrended VHF, UHF, and L-band strip charts from a pass recorded at Poker Flat on 29 May 1976. 


Courtesy: Radio Science, Volume 13, Number 1, pages 167-187, 
January-February 1978 AGU American Geophysical Union 


Electron Density Profile cintillation 
600 - t - + 800 . 
700 700 
600 600 
00 500 
E 400 = 400 
om & 
. 300 300) 
F-Region 2 ,__F-Regio 
Irregularities 20 200} 
E Layer —~>'™ , ey a 
0 es 0 : 
9 2 4 6 0 03 06 09g 
, [elicm>} x10" $4 index 


Results must be checked with real situations, from different sources. 

We chose LiveCQ as a source. 

Rene PE1L accepted to store all decodes from MAP65 spotters (all 2m band) in a file. 
We made an Excel sheet, with data sorted by date, spotter and spotted. 


Example: 18/08/2012 — DGOOPK — PE1L, data, pol and level graphs 
Note: MAP65 rotation is sum of spatial offset and up going and return Faraday rotation 


2 |date {rJute ({r)sig {*)pol {Jcall _[+spotter(~} 
18082012, 802, —-23—=—«SDIDGOOPK —_—PEIL 

18082012) 804) -25| {28|DG0OPK __|PE 
7) 18082012) BUR) -23) 448\/DGOOPK PE 
19082012, 943. ~-21'~~=«*A73|DGOOPK PEt 
18082012) 945 | -21| 8/DGOOPK _|PE 
qe0e2012,—953|+~=S=C«-28)~~=SCA DGOOPK PET 
f8082012| 9851 ——-22/—~—=BT|DGOOPK (PE 
18082012) 957 -23| \74\DGOOPK _|PE1 
18082012, 969-21 5BIDGDOPK _PEt 
18082012) 1001) -22| 58/DGOOPK _|PE 
18062012; 1016) _—-24/ 82/DG00PK PEt 
18082012) 1018) -20) 83/DGDORK —_PE’ 
18082012) 1020) -21| 86/DGOOPK _|PE1 
yeoe2012,1022/~=«=-20,~=~=S BS DGDOPK PET 
180820121024) -25| 117|DGOOPK —_|PE1 
18082012| 1026, ———-20,——«B7|DGNOPK PE! 
48082012) 1100) -25) 434|DGOOPK —_|PE1 
18082012) 1102, —~-22/~=S«13B,DGOOPK PE 
48082012) 1104) -23| 4146|DGOOPK —_|PE 
18082012) 1112) -23} 153|DGOOPK |PE1 
18082012 | 114) -22| 161|DGOOPK PE 
18082012) 1330) -24| 66/DGOORK _|PEt 
teog2012, 1508, 23) {1|DGOORK PE! 
18082012) 1510 -24| A7/DGOOPK _|PE1 
te0e2012) 1612; 25) 17IDGOORK PEI 
18082012) 1514) -23| 15|DGOOPK _ PE 
18082012) 1520) -22| B|DGDORK PEt 
18062012/ 1522, —-22/ 4|DGDOPK PE 
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lonospheric Absorption vs Frequency 
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Negligible on the 


higher bands and in 
night conditions. 
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Source:Radio Wave Propagation by Lucien Boithias, published by North Oxford Academic 


In 2 m JT65B decodes we see fluctuation of the levels, showing both 
medium term (4'-8’) ripple (2-3 dB) and long term (1-2 h) undulations (4-5 dB). 
ve APPILIPISIOPILEPELLS 
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Cannot be attributed to variation of attenuation. Most logical explanation 
is focusing or defocusing in curved layers of ionospheric waves. 


In 2 m JT65B decodes we see fluctuation of the levels, showing both 
medium term (4'-8’) ripple (2-3 dB) and long term (1-2 h) undulations (4-5 dB). 


Moon is wide 0,5 deg 


Our antenna beam 


is wide many degrees 
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Cannot be attributed to variation of attenuation. Most logical explanation 
is focusing or defocusing in curved layers of ionospheric waves. 


¢ In 1845 Faraday discovered that the plane of polarization of linearly 
polarized light, traversing a medium, can be rotated by the application of an 
external magnetic field aligned in the direction in which the light is moving. 


An electromagnetic wave, crossing the ionosphere, will rotate by: 
® = k* B* TEC/ f? (rad), with: 

B = Geomagnetic field component in Moon’s direction 

TEC = Total Electron Content of the path 

f = wave frequency 


Band dependence, with same B and TEC: 


50 MHz 90° 360° 2.25 turns 
144 MHz 10° 40° 90° 
432 MHz 191 479 10° 


1296 MHz 051 0S5 191 
Evidently, Faraday is a concern mainly in VHF 


Microwavers are concerned only by Spatial Offset 

Polar polarization is the angle between an antenna and earth's polar 
axis. 

Spatial offset between two stations is simply the difference 
between the polar polarizations of the two stations. 


For solving the algorithm we need sources for B and TEC 


Longitude: degrees 


From the web site of the British | a cae 


Geological Survey, introducing = 


Date; [2013-01-20 | (YYYY.YY OR YYYY-MM-DD) 


Lat&Long of station, 
Median Height of the ionosphere, | 
and Date, one obtains: 


Total field F (nTesla 
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We need B, Geomagnetic field component in Moon’s direction. 


Vector F is defined by -—Inclination and Declination. 
Vector Moon’s direction is defined by Azimuth and Elevation. 
For projecting Field F on the Moon's direction we need the angle FM 


between these two vectors. Formula: 


cosFM=cosl*cosD*cosEL*cosAz+cosl*sinD*cosEl*sinAz-sinl*sinEl 


B =F * cosFM 


TEC (Total Electron Content) is measured in TECU (TEC Units) = 10'® electrons/m? 


The number of TECUs represent the total number of electrons present in a cylinder of 
1 m* of section, crossing the ionosphere in the wave’s direction. 


We used data from the Royal Observatory of Belgium (ROB), in Dourbes, which 
publishes VTEC histograms with values every 15’, and archives each day of the year. 


The ionosphere cannot be defined by a number, Hi 
since its density varies with altitude. 
A useful schematization is representing it by a slab — , 
of uniform density. a ae 
This slab represent the transformation of the real ; wae 
ionosphere in an equivalent ionosphere . 


With two numbers we can represent an 


equivalent ionosphere. a ote t = -- _ Team 
The ROB (Dourbes) site gives both VTEC {1 ft \Seeep 


and Slab Thickness = ——_ SS SsSCFS‘—Fs 
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Global VTEC Map 


Total Electron Content (TEC) 2013-01-16 13:00:00 UT 
¢ TEC Longitudinal variation: Global trend quite 
regular and correlated to the local solar time 


¢ TEC Latitudinal variation e . 
The TEC value, varies non-linearly from the “= RC, = 


poles to the equator (geomagnetic) — 
With the algorithm representing this curve, Institute of Communication and Navigation, German Aerospace Center (DLR) 
introducing the Mag. Lat. of the station, on or red =|. TTT yd 

we find the correction of Dourbes VTEC. 


Slant TEC = See 
Crossing the slab obliquely there are more electrons. WacNeTiC LATTUDE /Nfaalabfenlsal se [ ecemere 
Instead of Vertical TEC we must use Slant TEC. TECU variation = 0,02*LAT2-2.5*LAT+95 
Path 
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TEC = STEC = Ka*VTEC 


With Earth radius=6367 km, lonosphere beginning at 
100 km height, and h=Slab Thickness 

Ka =(SQR((6467+h)?-(6367*cosEl)?)-SQR(64672- 
(6367*cosEl)?))/h 


ie 


¢ We now have the data for the complete formula. 
For 144 MHz, k/f*=1,14 with F in Gauss. 


¢ Wave plane rotation is controlled by these variables: 
LJ Angle FM between Geomag. field and Moon direction 


N hemisphere: 
cosFM ranges from 0 to -1 
S hemisphere: 
cosFM ranges from 0 to 1 


LJ TEC (constant or changing slowly, 100% to 30%) 


LI Moon elevation (slant passage Ka from 3.7 towards 1) 


We made an Excel sheet, and we got good congruence in the 
majority of cases analyzed. Example: 


[=e Declination(}| caic.Euiline 


. a 
29.994 
rae 


Ol 


1 TE 
Pan33wT | aa | oT | 9-0 41456 | 


| UTC | Localtime | Decimal [az() [El sretteaca| ¥TEC-ORBS | ¥TEC creewp 
24,38 a 
poss [ose [28.38 Ei) A 


[een tauae | Dectinationt 
0.953 
Fem [oe [Po | 


a a ee “41456 508,63 85,27 
Ee “DAT456 482,65 82,70 


Resumes | [| 


Having now confidence in the basic correctness of formula and correction 
coefficients, we proceeded to build a new Excel sheet, 

covering the entire common-moon period. 

Partial checks were possible using the LiveCQ decoded periods. 


Example: SP4MPB spotted by PASFPQ, total pass: 


| 2063] 50.65] 3] zo] osagse] parr] 4s] 
eel [ae | PD 


nS 
| hikm) | Ka | vTEC Orbs | Con —[¥TEC toc] STEC | cosFM | _Offset PT 
2 7) a 2 Tr) 


, 

it 
ai 1D a 
| co on |= | 00 |e } 
SPO; 


JOULES | 


| 


l 


P23] 36a] es] tz4o] 31.36] -0,6042] 
| aso] tea] a] go] 2,74] -0.6435] 
[tof] 12,82 


23,73 -O,6716 


t<] 
> 


21,40 
[20.60] 
a) 


\= 
5 

1 

JE 

| 


ih 


= 


STEC 


0, 

i [20 

a a es] 33 rm 
18.30 18.35 ; 63] 0.20 i 


222 
3 


3 f= J 


-0,3655 


(aboulate cl Pod. 


<= Decoded pol. from LiveCQ 


¢ SP4MPB was active from 13.58 to 14.42 utc (near sunset) 
¢ Inthis phase, TEC had a quick decrease. 


¢ Followed by a brief increase pre sunset, then decreasing 
from sunset to night. 


¢ Calculated and real trend are coherent. 
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¢ Spotter IK1 UWL (Band 144 MHz - Dec 19, 2012 — Moon 11.00 — 23.00 UTC) 
¢ All graphs computed for stations in a rose of directions 
qb=k*(F*cosFM)*(VTEC*corr*Ka)/f2 
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¢ Spotter IK1 UWL (Band 144 MHz - Dec 19, 2012 — Moon 11.00 — 23.00 UTC) 
¢ All graphs computed for stations in a rose of directions 
qb=k*(F*cosFM)*(VTEC*corr*Ka)/f2 
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Pol. trends as function of direction 


¢ Spotter IK1 UWL, (Band 144 MHz - Dec 19, 2012 — Moon 11.00 — 23.00 UTC) 
¢ All graphs computed for stations in a rose of directions 
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qb=k*(F*cosFM)*(VTEC*corr*Ka)/f2 
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QSB of JT65 decodes: 


Is caused by focusing or defocusing of our beam going through the waves of 
the windy ionosphere. 


Faraday rotation: 

There are three phases in a Moon pass: 

1 - In the first hours after Moon rise the rate of change of polarization is high. 
Causes: 

a) — change of angle FM between Moon direction and magnetic field 

b) — change in length of ionospheric crossing (slant coeff. Ka) 

2 —I|n the central part of Moon pass changes in angle FM and coeff. Ka 
balance each other, so polarization changes depend mainly from ionospheric 
evolution (of Total Electron Content) 

3 — In the last hours before Moon set the rate of change of polarization is high 
for the same causes of phase 1 


- Aspects of Weather and Space Weather in the Earth’s Upper Atmosphere: The Role 
of Internal Atmospheric Waves by Michael C. Kelly. 


- INGV istituto nazionale di Geofisica e Vulcanologia. 


- TOTAL ELECTRON CONTENT STUDIES OF THE IONOSPHERE John A. 
Klobuchar,, e t al Air Force Cambridge Research Laboratories L. G. Hanscom 
Field, Massachusetts. 


Te Potential of Broadband L-Band SAR Systems for Small Scale lonospheric TEC 
apping 

(Remote Sensing Technology Institute, German Aerospace Center (DLR) 
Oberpfaffenhofen, D — 82234 Wessling, Germany 


- Institute of Communication and Navigation, German Aerospace Center) 


- GEOMAGNETISM TUTORIAL Whitham D. Reeve Reeve Observatory Anchorage, 
Alaska USA 


- Frederick University, 7 Y. Frederickou St., Palouriotisa, Nicosia 1036, Cyprus 


- Electron density measurements of the plasmasphere — experimental observations 
and modeling studies 


- Cooperative Research Centre for Satellite Systems Department of Physics, La 
Trobe University Bundoora, Australia 


- Propagation Factors In Space Communications ( NATO) 


- Seasonal variations of storm-time TEC at European middle latitudes Royal 
Meteorological Institute (RMI), Belgium 


- Radio Wave Propagation by Lucien Boithias, published by North Oxford Academic 


By Giorgio Marchi IK1UWL and Flavio Egano IK8XTV - All rights reserved (August 26 - 2014) 


